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We have investigated the interrelation between the spin glasses and the structural 
glasses. Spin glasses in this case are random magnets without reflection symmetry (e.g. 
p - spin interaction spin glasses and Potts glasses) which contain quenched disorder, 
whereas the structural glasses are here exemplified by the homopolymeric globule, which 
can be viewed as a liquid of connected molecules on nano scales. It is argued that the 
homopolymeric globule problem can be mapped onto a disorder field theoretical model 
whose effective Hamiltonian resembles the corresponding one for the spin glass model. 
In this sense the disorder in the globule is self - generated (in contrast to spin glasses) 
and can be related with competitive interactions (virial coefficients of different signs) and 
the chain connectivity. The work is aimed at giving a quantitative description of this 
analogy. We have investigated the phase diagram of the homopolymeric globule where 
the transition line from the liquid to glassy globule is treated in terms of the replica 
symmetry breaking paradigm. The configurational entropy temperature dependence is 
also discussed. 



I. INTRODUCTION 

The existing dynamical theory of the glassforming overcooled liquids and polymers are mainly applica- 
ble to the relatively high temperature area. The well-known mode - coupling theory (MCT) predicts a 
dynamical phase transition at the critical temperature Tc where the overcooled liquid becomes nonergodic 
and stays in this state upon the further cooling. It is significant that according to MCT, the characteristic 
times of the a - relaxation or the viscosity exhibit a power low divergence: Ta{T) oc ri{T) (x \T — Tc\ ai 
a critical temperature T = T^. It has been realized now that such types of singularities are generic for a 
broad class of mean - field (MF) models and its appearance is an indication for the role of the activation 
processes which are neglected in MF - models ||^ . This means physically that the topography of the free 
energy landscape in the space of the coarse grained variables should be taken into consideration. At low 
temperatures the free energy landscape determines more or less the dynamics of the system. The size of 
the free energy barriers between the metastable states, however, determinates the rate of any activation 
processes at the low temperature regime T < T^. 

At the so-called Kauzmann temperature Tk, which is usually Tk < Tc a genuine thermodynamic glass 
transition or so called Gibbs - DiMarzio transition is expected to take place ||^ . At the Gibbs - DiMarzio 
(or Kauzmann) temperature Tk, the system is preferably in a thermodynamic unique configuration since 
the height of the free energy barriers grow with the system size. 

The earliest analytical approach to the problem which captured these two aforementioned aspects has 
been worked out in several papers by Kirkpatrick, Thirumalai and Wolynes in the late 80 's 0. In these 
papers which study the phenomena by making use of the p - spin interaction spin glasses (at p > 2) and 
Potts glasses with more then four components, the following conclusions have been drawn: 

• The predicted phase transition temperature Ta obtained by the dynamical theory (actually equiv- 
alent to Tc of the MCT based on the Langevin dynamics in the mean field limit) is higher than the 
Kauzmann Tk obtained by the static theory (or by the ordinary replica theory) . 

• As the temperature decreases (starting from the paramagnetic phase) the metastable minima appear 
first exactly at Ta. 



In the intermediate temperature regime Tk < T < Ta many metastable states are separated by high 
barriers. Therefore these states may have a long lifetime. The activated transport is the typical 
process in this temperature range. 
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• The overlap order parameter within the one-step rephca symmetry breaking (1-RSB) scenario un- 
dergoes a discontinuous jump at Tk- This reminds to a first order phase transition even so ther- 
modynamically (e.g., in the specific heat) the transition seems to be of second order. The authors 
called this class of phase transitions random first order phase transition. Crisanti and Sommers 
confirmed essentially the same type of behavior for the p - spin spherical model Q , which is dif- 
ferent from the Sherrington - Kirkpatrick model, and is shared by the other spin - glass models 
without reflection symmetry 

Indeed it has been argued often, that many connections exist between the behavior of the structural 
glasses and spin glasses without reflection symmetry (|,^,^ . Most of these properties mentioned above for 
the spin glasses can be found in the context of structural glasses too. However, the important difference 
between the spin and structural glasses is that the spin glasses models contain a quenched disorder 
already in the Hamiltonian. In structural glasses the Hamiltonian is a regular function of the particles 
coordinates and does not contain disorder. Nevertheless the free energy in both systems resemble each 
other and posses similar properties. Thus, and disorder is in a sense self - generated and develops itself 
during the cooling (or glass transition) process. 

The properties and possibilities of self - generated (or self -induced) disorder have been already discussed 
in the framework of rather special spin - models These models explicitly involve pseudo - random 

numbers, i.e., the spin -spin couplings, although deterministic, oscillate very rapid and can be considered 
as a (pseudo) - random variables. These models, however, provide some spin analogies of structural 
glasses but still are quite different from structural glasses [|| . 

Furthermore, the ideas and the methods of the mean - field spin - glasses ^|-|l0| have been successfully 
used to study the freezing states in a heteropolymer globule JTl-p^. The corresponding model considers 
the conformational behavior of a polymer chain with randomly quenched interactions between monomers. 
These random copolymers posses always a collapsed (or globular state) depending on temperature and 
strength of the randomness. It was successfully shown, that the disorder results in a globular and 
eventually in a glassy globular state. The corresponding freezing is characterized by a transition between 
two phases: one phase is characterized by many accessible configurations, while the other is dominated 
by only a few of them, i.e., an eventually frozen state. In this context the heteropolymer freezing usually 
serves as a simple "toy model" for the protein folding phenomenon [ [T^ . This relationship between the 
freezing and folding phenomena is caused mainly by the fact that in both cases only a single conformation 
(which is in the context of protein folding called the native state) dominates. 

In the following we are going to investigate a similar problem in homopolymer globules in poor solvent. 
It is well known that homopolymers in poor solvent form globules due to an attractive second virial 
coefficient which are then stabilized by the repulsive third virial |]l5| . The relevant parameter for the 
globule formation is the second virial coefficient which is measured by the relative distance from the 
so-called Q - temperature, where the homopolymer takes Gaussian conformations. Close below Q the 
polymer collapses and the deeper the temperature the more dense becomes the globule. The intuitive 
limit is then a dense liquid globule. We will show nevertheless, that even at lower temperatures and 
beyond a certain density the globule freezes and forms a glass. Homopolymer globules show a structural 
glass transition. Thus they will form a "nano glass" . 

It is most interesting to note that in recent Monte - Carlo simulations of homopolymer globules, Milchev 
and Binder JTot have found a dramatic decrease of the acceptance rate of the moves which in its turn 
suggests to the growth of the characteristic times. They have seen also a pronounced density fluctuations 
which appear in the center of the globule and spread over the entire globule as a temperature is lowered. 
This results have been interpreted in terms of a glassy type of transition. These investigations have been 
restudied in papers by Kreitmeier et. al p^-p^ where a similar dynamical behavior of the homopolymer 
globule has been verified once more and was generalized for the cyclic deformation regime. Moreover, a 
glass transition in a homopolymer globule (for the chain length < 27 on the 3x3x3 cubic lattice) 
has been discussed by computer simulations |2^. These authors found that the conformation space of 
the model consist from one huge valley and a big number of small "chambers" . These disjoint domains of 
the conformational space are responsible for the ergodicity breaking and freezing transition. This might 
be important for the final processes of folding dynamics in the protein molecules because the speed and 
reliability of the folding depends on whether the native state belongs to the same conformational space 
domain or not 

The thermodynamic behavior of the homopolymer globules have been studied also by discontinuous 
molecular dynamics simulation |2iy22f| . By making use a simple off-lattice model for chain length = 64 
the authors showed the existence of the first - order liquid - solid - like transition. The transition occurs at 
temperatures below the coil - globule transition temperature and has been detected by the heat capacity 
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peak as well as by the Lindemann parameter jump. It is interesting that there is a qualitative similarity 
between homopolymers and proteins where the transition from the molten globule to the native state is 
also of the first order [^,|^. In ref. |^ a new "expanded ensemble" Monte Carlo algorithm has been 
introduced which helps to overcome the density slowing down in the globule state. For the chain length 
up to = 512 the authors have seen the bimodal distribution in the number of contacts per monomer. 
This is again a clear indication of the first - order liquid - to - solid transition of the collapsed globule. 

As already mentioned, most important is that the globule formed from ordinary homopolymers does not 
contain any quenched disorder. If this system forms glassy states the disorder has to be self - generated. 
Nevertheless, to form glassy states certain frustrations are necessary. These may have its origin in the 
interplay between attractive interaction (negative 2"*^ virial coefficient v), repulsion (positive 3"^ virial 
coefficient w) and the constraints imposed by the chain connectivity. These frustrations might be already 
sufficient to provide the structural glass transition from a liquid globule to a glassy globule. 

The main purpose of the present paper is to investigate these problems in more detail by analytic 
means. We will indeed show that the linear homopolymer in the condensed globular phase undergoes a 
genuine structural glass transition of the similar type as in the p - spin - interaction spin - glass with 
p > 2 or Potts glasses models. We stress once more, that the homopolymer is a pure system and does not 
contain quenched disorder. Nevertheless the n - component field theory formulation for linear polymers 
p6t provides already a natural and solid basis |^^. As we will show below, the homopolymer globule 
is a generic system where frustrations rather than quenched disorder trigger the glass formation. As a 
starting point we employ the field - theoretical description of the self - interacted polymeric chain in a 
poor solvent |^,^. So far, in ref. |^,^ it has been shown that the Laplace transformation (with respect 
to the chain length N) of the polymeric correlation function can be associated with the corresponding 
correlator within the n — > limit of the n - component '0a, V'a " field theory (where a = 1, . . .n). We 
argue instead that the same n - component field theory can be mapped onto some kind of disordered one 
- component model where the imposed random field is a non - Gaussian (a color noise) and its statistical 
moments are determinated only by virial coefficients v,w,. . ., etc. In this case the index a acts as a 
replica index for spin glasses |^J^ or random magnetics ||l0|| . This mapping has been suggested for spin 
glasses already three decades ago [ pOpl] ] but was never actually used. At this point we should add a 
general remark: Throughout the paper we use the large globule approximation, i.e., we consider the chain 
length N to be very large. This is necessary to avoid additional complications close to the surface of 
the globule. Intuitively it is obvious, that monomers close to the globule surface experience a different 
dynamic environment, such that different mechanisms of freezing processes take place. We will come 
back to this point in a later publication. 

II. THE HOMOPOLYMER GLOBULE AS A DISORDERED SYSTEM 

A. The field theory representation 

In this section we will provide the basic formulations and the field theory for polymers in poor solvent, 
i.e., with attractive interactions. To do so, we start from usual continuous description of a homopolymer 
chain of the length N in the poor solvent. Its Hamiltonian has the following form 

N 2 ^ 

"^'^'^^^^ I ^'(^) -^41 dsds'S{ris)^r{s')) 



N 

+ 11 J dsds'ds"S (r(s) - r(s')) S (r(s') - r(s")) , (2.1) 



where we allowed already that the 2"^^ virial coefficient is negative, v < and the 3''^ virial coefficient 
w > 0. Usually the second virial coefficient is measured by its distance from the & temperature, i.e., 
\v\ ~ a^|T — Ql/Q = a^T. The third virial coefficient is always of order of w ~ 0{a^). The density of 
the globule is easily estimated by p ~ which provides a naive limit of the temperature at pa'^ = 1. 

Indeed at sufficient low temperatures the globule is dense, and the entropy term (first term) becomes less 
and less important. The remaining entropic fiuctuations are only important on length scales ^ ~ a/r. 
This observation will allow us later the safe use of corresponding mean field theories. 

The next step is |2^j2^,^ (see also Appendix A) to employ the the Laplace transformation of the 
polymeric correlation function S(ri, r2] N) 
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S(ri,r2;Ai)= J diV S(ri, ra; A^) exp(-^7V) 



which can be associated with the n hmit of the -0^, tp^ field theory as foUows: 

Y[ I?VaV'i(i"i)V'i(i"2)exp|-iJ„[7/^;^]| , 



(2.2) 



(2.3) 



where ip — {ipi, 1(^2, ■■■ , 4'n\ is a n - component field, /i , the chemical potential, which is conjugated to 
N and the replicated Hamiltonian is 



1 /• " 



2d 



w 
3!8 



H 3 



(2.4) 



As usual the vector field ^p corresponds to the polymer density in the usual manner, i.e., p cx (ipl) The 
relationship between p and N has the form 



TV : 



J d'^rid'^r2E{ri,r2;p) 



J d<^rid<^r2^{Yi,V2;p) 



(2.5) 



The attractive interaction term of the order -0* changes the behavior of the field theory. Correlations of 
the self avoidance are no longer important and we must take care on the balance between the attractive 
and repulsive forces. Although we mentioned above the relative unimportance of the connectivity term, 
we have to take track on it as well. The careful analysis below shows that it provides at sufficient 
low second virial coefficients significant contributions on small scales (^ ~ C'(<^)) contributions only by 
connectivity. These are in part responsible for frustration. 



B. Mapping onto a random model 



Now we are going to map this field theory onto a random system. It can be shown (see Appendix B) 
that the free energy of the globule state Fq\ can be interpreted as the free - energy of a one - component 
random model with Hamiltonian 



where the random field t(r) is non - Gaussian with its generating functional of the form 



(2.6) 



${p(r)} = expj-y" 



exp 



dd-rt{r)p{r) 
d'^rp^ir) - 



w 
3!8 



d'^rp^ir) 



(2.7) 



and where the bar means the averaging over t{r). It is interesting to underline that only for this combi- 
nation of signs of the virial coefficients {v < 0,w > 0) the even central moments t{ri)t{r2) ■ ■ ■ t{r2m) (see 
eqs.(|B^) - (Bll)) are positive, as it should be for the real field t(r). The aforementioned mapping takes 
the form 



Fgi 

where the replicated partition function is 



lim — log Z", 



D^pe:cp{-Hn} 



(2.8) 



(2.9) 



In the present paper we are going to use this analogy between the homopolymer model (which is a "pure" 
model, i.e. does not include a quenched disorder in its Hamiltonian) and a random model (p.6[). 
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C. Legendre transformation 



As a next step we should go to the two - rephca - variables, Qabi^), or Parisi overlaps |^ for polymers 
p7| . One can merely implement it using the Legendre transformation of the interaction part of the 
Hamiltonian (2.4) which can be represented in the form 



b^bc^ca 



a,h=l 
Z 



" >— 1 ■ a,b=l 

« n 

■ j d'^r UabUbcUcdUda + 0{uli^), 



(2.10) 



where the pair field Uob(?') = '4>a{T)'4'b{f') and where we have also kept the 4"^ virial coefficient z. Let us 
introduce the integral transformation 



e^]i{K[uab]} 



J n DQ^dir) exp J W[Qab] +J d^rj^ Qab{r)uab{r) \ , (2.11) 



and find W^[Qa b] in the form of a functional expansion. For this purpose one should use the saddle point 
method in eq. ( ^.llD which can be carried out in the same spirit as in references [p3^,p4[. This results in 
the Legendre transformation with respect to the extremum field Q^jj 



/n 
d'r ^ Qabir)uab{r). 

a,b=l 



(2.12) 



As a result [B5| we get: 



OUab 



(2.13) 



6W 
SQcdir) 



(2.14) 



By making use of the expansion ( 2.1C ) in eqs.( ^.13 ) and ( 2.14 ) one obtain for ^^[Qab] the following 
expression: 



W[Qo.f,ir)] 



d''rTr(Q2) 



4w 
2 



w 

M 



d'^rTr(Q4)+o(g5 



After the transformation (2.11) the replicated partition function (|2.9|) takes the following form 



r) exp 



Trlog \ [^i - ^'^^]Sab - 2Qabir) }> + W[Qabir)] } , 



(2.15) 



(2.16) 



So far only mathematical identities have been used. Nevertheless the use of the overlap variables Qab 
allow to detect completely different correlations as the classical 0{n) field theory in the limit n ^ for 
self avoiding walks. They will allow to probe for a more complicated phase space and provide information 
on the presence of glassy type correlations in the globules. Wit h th is in mind it appears instructive to 
express the generalization of the polymeric correlation function (2.3) 

Sab(ri,r2;Ai) = {ipa{ri)A{r2)) (2.17) 

in terms of overlaps Qabif)- For this end we add the source-term '4'a{T)ha{r) in eq. (^.ll ). After substi- 
tution in (|2.9|) and integration over -0 one gets 
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^ah(r,r';^) 



1-Gn 



2Q 



(2.18) 



where Gq^ = IJ-— {a? /2 d)V\ The corresponding polymer correlator ( ^.3| ) is nothing but the Sii(r, r';/i) 
element of the matrix ( ^.18| ) . 

The correlator Sah(r, r';/x) measures the probability to fin d cha in configurations started at r in the 
replica a provided that it ends at r' in the replica b. From eq.( 2.18| ) this correlator also can be seen as a 
scattering amplitude of a free "particle" with the Green function Go on the "scatterers" whose density 
is described by Qcd- 

The present representation ( 2.16| ) is very promising and reminds on the corresponding expr essio ns for 
the spin glass models IHj^,^,^. The striking difference between these two cases is that in eq.( 2.16 ) only 
the pure model parameters (Kuhn's segment length and virial coefficients) are involved. As a result the 
representation can provide a good starting point for the phenomenon of self - generated disorder which, as 
we believe, is behind the structural glass formation. 



III. MEAN FIELD TREATMENT 



A. Landau-type expansion 



In order to simp lify the mean field treatment of the integral ( p.l6| ) let us expand the effective Hamil- 
tonian in eg. ( 2. 16 ) to the 4"^ order. We then obtain by this procedure 

n DQab{r)ex:pl- ^ d''r,d''r2T'^^\ri,r2)Qab{ri)Qba{r2) 

a, 6—1 k a, 6—1 

J d'^rid'^r2dSr(3Hri,r2,r3)gafc(ri)Qfcc(r2)Qca(r3) 

a,fe,c— 1 

/rfSd''r2d'^r3d''r4rW(ri,r2,r3,r4)Qab(ri)Q6c(r2)Qcd(r3)Q<ia(r4)l, 



(3.1) 



a.b,c,d—l 

where the coefficients 

r(^Hri,r2) 

r(3)(ri,r2,r3) 

rW(ri,r2,r3,r4) 



-6{ri - ra) - Go(ri - r2)Go(r2 - ri) 



r6{ri - r2)5(r2 - r3) - Go(ri - r2)Go(r2 - r3)Go(r3 - ri) 



w 

M 



S{ri - r2)(5(r2 - r3)S{r3 - r^) 



Go(ri - r2)Go(r2 - r3)Go(r3 - r4)Go(r4 - ri 



(3.2) 



As often in the MF-theories the order parameter Qab does not depend on the spatial coordinate r. As it 
is custom in the MF-theory of the spin glass models [^-p^ , we decompose the Parisi matrix Qab in the 
following form 



lab = (g - f)Sab + f + Aa&. 



(3.3) 



In eq.( |3.3| ) the symmetric part Rab = (? ^ f)Sab + .f, with the diagonal q and off-diagonal / elements, 
describes the replica symmetrical (RS) solution ]^-|l0|]. The matrix A^fc equal zero for a — b and is 



responsible for the replica symmetry breaking (RSB) |^|-|1C|. The use of the decomposition (3.3) in 
eq.(3.1) allows in the MF- approximation to separate the total free energy into RS - and RSB - parts. 
During the calculation of the traces in eq.(3.1) it is convenient to use Parisi's representation of A^fc by a 
function A(a;), where < j: < 1. Then for the free energy is found to be 



lim 4- F{Qab} = /Rs{g, /} + /rsb{9, /; A(a;)}, 
riV 



(3.4) 



where RS-free energy 
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/Rslg, /} = A{q^ - f) + B{q' ~ 3qf + 2f) + C{q* - 6q\f + 8qf - 3/^) 
and the RSB -free energy 

r-i r /•! 



/RSB{g,/; A(a;)} = -wi /" dxA\x)-W2 f dxA{ 



(x) 



-W3 / dx 
Jo 



xA''{x)+3A{x) I dyA\y) 




W4 



dxA{x) 





+W5-I4 / dxA{x) ( dyyA^{y)+ f dxx'^A^{x) 
Jo Jo Jo 



+2 / dxA'^ix) 



dyA{y) 



1 -|2 

dyA{y) 



2 / dyA{y) [ dzA{z) 

Jo Jy 



The coefficients m eqs.(|3.5[) and (3.6) have the following form (see Appendix C for the details) 



\v\ V 27ra 



4-d 

'"^ r (g^) 

2^— 



C 



2 z 



27ra^ 



d/2 -p /'8^\ 



(3.5) 



(3.6) 

(3.7) 
(3.8) 
(3.9) 



wi = A + iB{q- f) + QC{q- ff 
W2 = -35/ - 4C/'' 
W3 - -S - 4C(g - /) 
W4 = -4C/ 
W5 = -C. 



(3.10) 
(3.11) 
(3.12) 
(3.13) 
(3.14) 



The minimization of /rs leads to the RS-solution, q^ and /„, whereas the maximization (as it is the 
case for the spin glasses [p|-p^) of /rsb results in the RSB-solution in terms of the overlap matrix A(a:). 
The coefficients in the Landau - expansion of /rsb depend also from the RS-solution. 

In order to take into account the spatial correlation in RS - sector (see Sec.V B), it is convenient to 
assume that the variable q and the coefficients A, B and C are weakly k- dependent. In this case the 
Landau - expansion is more involved and is given in Appendix C. 



B. Density in terms of Qab 

The essential issue is to express the globule density p in terms of the order parameter Qab- This 
allows to detect any glassy features inside the globule and to distinguish between the liquid and glassy 
phase. Moreover it will show any unusual properties of the phase space. As mentioned already in the 
introductory remarks the density is determined by the '0 - fields. Especially in the MF-approximation 
the local monomer density is given in terms of the grand canonical partition function as 



N 



P(r) 



ds 



J d'*rS(r;s)S(r-r';7V- s) 
JddrE{r;N) ' 



which should be supplemented by the normalization condition (mass conservation), i.e., 



d'^rpir) = N. 

By making use of eqs.( p.l5| ) , ( 3.16| ) and after Laplace transformation we find 



(3.15) 



(3.16) 
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(3.17) 



where S(r; /i) is the Laplace transformation of ^(r; N) ( see eq.(P^)). Then the equation which determines 
the chain length becomes 



and takes the form 



N 



N = 



d_ 



logjy d''rS(r; ^) 
/ d'^r J d'^r'Eir'; ^)S(r - r'; ^) 



/ d''rS(r; N) 
In the MF- approximation eq.(^.3|) reads 

S(r) = <(rX(0), 



(3.18) 



(3.19) 



(3.20) 



where ip™^{r) is the MF-solution for ■ipi{r). Combining eqs. ( |3.20| ) with (|3.19[) and (^.16|) leads to the 
expected result 



p(r) = [Vr(r)] 



(3.21) 



In order to express finally the densi ty in ter ms of Qab we recall that the pair field u5"/= [V'5°''(i-]Il( see 
Sec.IIC). Combining this with eqs. ( ^.14 ) , ( 2.15 ) and taking into account the decomposition (3^) one 
obtain for the density the following expansion 



4w 



8 /w 



\v\ 



4w 



dxA'^{x) 



V 



8 /'w^ 



V 



3q / dx^\x)+ / dx (xA^ix) +3A{x) / dyA^iy) 



(3.22) 



In eq.(3.22) we have used Parisi's representation of Aab and taken into account that the off-diagonal 
element vanishes, / = (see below). The "singularity" \v\ ^ in eq.( 3.22 ) is spurious, as we will see in 
the next subsection in the MF-approximation the value of q becomes \v\ dependent itself, i.e., q cx |up 
and A{x) oc jvp, so that p ^ at w ^ as it should be close to the 8 - temperature. 



IV. THETA POINT REGIME: COIL-GLOBULE TRANSITION 



The question which must be resolved first is the ordinary coil globule transition. In any case the 
pre sent general approach should reproduce the physical properties of the standard coil globule transition 
|p6| , ^ . To do so we investigate the system close below to the Q - temperature T < Q . In ref. ||28| , p9| 
the standard 0{n) - field theoretic formulation (see Sec.IIA) has been used in order to treat this problem 
beyond scaling. The method of pair fields, which is a simpler version of the present formulation p7|| has 
been developed and applied for the coil-globule transition. At that earlier paper of one of the present 
authors the third virial coefhcient has not been taken into account, which corresponds to an expansion 
around the Q point regime. Here, the Legendre transformation method from Sec.IIC allows easily to take 
into account an arbitrary number of virial coefficients. 

In the theta point region (i.e. at T < 8) the globule conformations are very close to Gaussian form 
so that for the chemical potential one can expect the following scaling: /i — fio/N. In this case it is 
convenient to rescale the virial coefficients \v\,w and z in the following way 

\v\ 4-d 

iN^, (4.1) 



a' 

r3-d 



y = (4-2) 



Z 8-3d 

t=Z^N^^ (4-3) 
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which shows the upper critical dimensions of the different terms in the virial expansion. After that the 
scahng form of the coefficients (3.7) , (3.8) and 



A- 



B 



C 



N— 



AN— 



-) 

2tt J 



is given by 



3 a 



2d 



N- 



^3d 



d 

2 t 

3 a? 



4-d 

Mo' - 










6-d 

2Mo' J 




, ^ y/2r(8^)- 


y2TT J 




3a*o' - 



(4.4) 
(4.5) 
(4.6) 



In the theta point regime \v\ = 0*^(1 — T/Q) and N >> 1, so that atd^S xwl, and 
t << 1, i.e. the forth virial coefficient becomes irrelevant. 

In the present regime only RS-solution makes physical sense, since n o ot her solution than the onset of 
the liquid globule can be expected. Thus we minimize RS-free energy (3^) with respect to q and /. The 
resulting solution reads 



/,„ = 



-3B + y/{3B)^ +32\A\C 



8C 



so that the RS-free energy becomes 

fns{q}^Aq^ + Bq^ + Cq\ 



(4.7) 
(4.8) 

(4.9) 



Let us consider the possible second order phase transition and impose the following conditions: A < 
0,B > and C > 0. In the vicinity of the transition point the coefficient A becomes small \ A\ << 1, and 
the order parameter takes the value 



3B ' 



(4.10) 



As is see n from eqs. (4.4) and (4.5), in d = 3 the order parameter scales as qm oc 1/N. Thus, as it follows 
from eq.( 3.22| ) the density scales as p = 4gm/|'i;|, which means physically p = l/N^^'^ (note also x « 1 
and \v\ cc 

Therefore we reproduce the correct scaling for the density, which is found also from 

naive scaling. 

We can also obtain the transition line in the \v\ - w plane, which will be the first step towards a 
more general phase diagram spanned by the virial coefficients. This line is defined by the conditions 
A = 0, S > and C > 0, which again for three dimensions = 3 yields 



Mo 



y_ 


> 




^3/2 


r(i) 










yl 


> 




^3/2 


r(l) 








6Mf 



(4.11) 
(4.12) 
(4.13) 



It is interesting that if eq.(4.11) valid then conditions (4.12) and ( [f.l3|) merge and conve rt int o the 
globule stability condition: y > 167r^/27. We can eliminate fiQ in eq. (|4.11| ) by combining eq.( 3.18 ) with 
the polymer correlation function 



a±i.2 
6 



p. -2q„ 



(4.14) 



We recall that in eq. (4.14) /i = fi^/N and q^ = q^-^/N . The result of the combining reads 
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MO = 1 + 2^-2^50^ . (4.15) 



By making use of eqs. ( p^ ) and ( pHl ) in cq.( [1.15[ ) for the transition line, one gets 



where the globule stability condition y > 167r^/27 is implied. 

For completeness we check for the possibil ity of a first order phase transit ion. Th e necessary conditions 



for this are:A > 0, B < and C > (see eq.(O)). It is simple to see by eqs. (|4.4[ - ^!^ ) that these conditions 



are contradictory. This means that within our MF-approach only the second order coil - globule phase 
transition is possible which is in accordance with the well known result . Therefore the present field 
theoretic formulation is able to reproduce the standard coil globule transition as the replica symmetric 
solution at conditions close to the Q - temperature. 

V. DEEPER IN THE GLOBULE STATE: LIQUID VERSUS GLASSY REGIME 

A. RSB solution in the globule 

In the following section we investigate the possibility of replica broken solutions deeper in the globular 
state. The globule density from naive scaling is given by pa^ = a^\v\/w = r and has natural limit at r = 1. 
Physically this limit corresponds to a dense globular state without any solvent inside. At temperatures 
below the coil - globule phase transition, but still far above r = 1 the system is usually characterized 
by a monomer - monomer correlation length ^ < i?Gi oc N^^^. In this case the chemical potential and 
the density are no longer N- dependent. It can be seen from dimensional analysis and simple scaling 
arguments |2^,|2^,^ that the chemical potential scales as/ioc|wp/wocT^, the density pci \v\/w oc r 
and ^ (X l/\v\ cx r . 

In this regime fluctuations can be still important unless a < ^ << Rgi, where the MF solution, which 
we have discussed in Sec. Ill, becomes valid. Here one can expect that because of competitive interactions 
(negative second virial coefficient versus positive third virial one) and the constraints imposed by the chain 
connectivity only a few conformations are dominated. This could manifest itself as the glass transition 



long known for spin - glass models [|-10| and heteropolymcrs [Q. Formally speakin g, th is transition 



shows itself as a nontrivial solution which maximizes the RSB - free energy functional (3.6) 



Before doing this extrcmization let us find first the corresponding e quati on for the chemical potential 



In the RSB - case we should calculate the Sn element from equation ( p.l§| ) and substitute it into ( 3.1^ 
We use here the so called one step replica symmetry breaking (1-RSB) scenario, which is generic for the 
glass transition in p - spin model random - energy model js^ and random heteropolymcrs [Q. 

Then Parisi's function A(2;) is defined by only two parameters 

A{x) = ( ° ' ^ ^ (5.1) 

The interpretation of a and the break point xq is the following Q . Within the 1-RSB scenario all replicas 
are grouped into two clusters so that 1 — xq is the fraction of the replicas which belongs to the overlap 
cluster with the overlap strength a. The rest of the replicas does not overlap. After this simplification 



the inversion in eq.(2.18) can be done analytically (see eq.(AII7) of ref. ||39|]). After a straightforward 
calculation one gets 

2ii(k;u) = ^ ^ — — (5.2) 

xo [^k^ + ^l - 2q - 2(1 - xo)a] xo [§3^2 + ^^2(7- 2a] 



Insertion of eq. ( 5^ ) in eq. (|3.18| ) simply yields 



- 2q„, - 2(1 - xo)a„^ = ^ { ] : (5-3) 



where qm and am are the solutions which extremize the free energies (|3.5| ) and (|3.6|) correspondingly. 
With 1-RSB assumption (see eq.(|5.1|)) the free energy (|3.6D becomes 
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/rsb(o-, xo) = - wi(l - xo)a^ + |w3|(l - xo)(2 - xo)cr^ 
- |w5|(l - a;o)(3 - 3xo + xl)a^. 



It is convenient to represent the chemical potential in the form 

|2 



V 



si\v\,w) 



(5.4) 



(5.5) 



where s(|w|,w) is a function of \v\ and w. For convenience we switch to dimensionless variables (by 
keeping for simplicity the same notations): 



w 



w 



,9 



z 



a" a" a" 

,3 , A D„3 , D /-<„3 



Aa-" ^ A , Ba-" ^ B , Ca-" C 
We also introduce the reduced (with a bar) values: 

A^A\v\ , B = B\v\' , C = C\v\' , = ^ 



(5.6) 



Wl = Wl \v\ 



W3 = W3|W| 



W5 = W5|i;| 



After that, the equation (5.3) for s(|w|,w) (or for the chemical potential) takes the compact form 



- = 2g„ + 2(1 - XQ)ar, 
s 



where 



-3B+ -^(35)2 +32|A|C 
"8C 



(5.7) 



(5.8) 



(5.9) 



3|w3|(2 - xq) + V(3|w3|(2 - xoj)^ - 32wi|w5|(3 - 3xo + 4) 
8\w5\{i - 3x0 + xl) 



(5.10) 



Wl = ^ + 2Cq, 
W3 = ~B - ACq, 
W5 = -C 



m 7 



and the reduced coefficients reads 



A = 2- 



27r 



3/2 



,1/2 



5^1 

3 



(5.11) 
(5.12) 
(5.13) 



(5.14) 
(5.15) 

(5.16) 



Eqs.(5^) and ( 5.1C ) are the result of extremization of eqs.(49) and (5^) correspondingly. Throughout 
the remainder of the paper we will retain the conditions : A < 0, B < 0, C > 0, wi > 0, W3 < and 
W5 < 0. This assures that the nontrivial solution a^n shows up via a first order phase transition, in a 
similar manner as in p - spin spin glasses or the random - energy model |38|. On the coexistence 
line between liquid and glassy phases additionally fnsB{o'rm xq) — and this leads to the corresponding 
equation 



Wl|W5| 



(2-xo)2 



|w3|2 4(3 - 3X0 + a;2) 



(5.17) 



11 



The equation for th e re duced globule density, c = p\v\, can be easily obtained from eg. ( 3.22 ) under the 
1-RSB assumption (5.1). The calculation yields 



4w 



Aw 



™ /I 2\ I 3/1 \2 

— {l~Xo) + -a^^l-Xo) 



—a^{l~xo) + j^[ — + - 



(5.18) 



Eqs.(5.8) - (5.16) for the function s(|w|,w) can be solved numerically at given values of the forth virial 
coefficient z and break point xq. After substitution of this solution s{\v\, w) in eq.(5.17) we arrive at the 
equation for the coexistence line in the plane of \v\ and w. By changing xq one can obtain a whole set 
of xo - isolines. We will give the corresponding numerical solution in Sec.V C but before we analyze the 
validity of the MF - approximation which has been given above. 



B. Role of fluctuations in RS - sector 



As have been mentioned above, the MF-solution is valid when the fluctuations are negligible. Generally 
speaking, this should be required for RS - and RSB - sectors of the replica space. In the present paper 
it is not our intention to consider fluctuations in the RSB-sector, which is rather involved, and we leave 
it for a future publication. 

In the RS - sector spatial fluctuations are described by the correlation function (C13) (see Appendix 



C). It is easy to calculate from eq.( pl3 ) the radial distribution function g{r) = ATnP{Aq{r) Aq{0)) at 



eq.(C13), the radial distribution function reads 



3. Again we turn to the reduced variables (5.5) - (5.7). After the inverse Fourier transformation of 



q(r) — ^= exp l — 3|ti|'\/ = r 
y\ ^ 2X I ' ' V 



(5.19) 



where 



X = A + 3q,nB + 6ql,C 



and 



Y 



1 

24tt V 2 



3/2 



(5.20) 



(5.21) 



It is important to note here once more, that the quantities X and Y do not have any additional |ti| - 
dependence. 

We estimate now the Ginzburg parameter, ea, (see e.g. ||40| ]) as the ratio of g{r) at its maximum to 



9m ~ 9ml''^l ■ Then for the Ginzburg criterion we have the following equation 

1 



\v\' = 



(5.22) 



QeeGq?,^^J X{\v\,w)Y{\v\,w) 



where e is the Napier number. Eq.(5.22) for a reasonably small Sg represents a line in the (|w|, w) plane 
which separate fluctuating and MF regimes. From now on we will call these two regimes as the liquid 
globule and the glassy globule correspondingly. 
From eq.( 5.19 ) the correlation length is given by 




(5.23) 



which is qualitatively in line with the standard result |26 2^,Q and the correct scali ng. 

From the Gaussian approximation for the effective Hamiltonian in the RS case ( C10| ) it is easy to 
calculate the corresponding RS - free energy, which takes the following form 
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fRs{qm} = \V\' [Aql + B<t + eta] + | / % ^ X + ^Y] 



(5.24) 



The last integral in eq.(5.24) diverges at large k. This ultraviolet divergence is of no significance as soon 
as t he \v\ - dependence is the only one we are interested in. For the second derivative of the integral in 
eq. ( |3.24 ) with respect to |up one has 



'1^ 



T fx 



I [\v\^X + n^Y\ 167rVr 



-v\ 



(5.25) 



After that the expression for the RS - free energy takes the form 



-v\ 3/2 



(5.26) 



We ascribe this branch of the free energy to the liquid globule state and will consider it in more details 
in the next subsection. 



C. Numerical calculations: chemical potential, phase diagram and configurational entropy 



Finally we are going to compute the phase diagram for the polymer globule. We recall here that the 
first step carried out earlier in this paper, i.e., the RS - solution corresponds to the classical coil globule 
transition. Here we are now in the position t o ca lcu late f rom the 1-RSB free energy the transition to the 
glassy state. The numerical solution of eqs.( ^.8|) - ( 5.16 ) at z = 9 and xq = 0.88 is shown in Fig.l . As 
can be seen the function s(|u|,w) depends linearly from w and almost does not depend from \v\. This is 



in agreement with well known result [E9U39]: ^ cx |u| /w (see eq.(5.5)). We have a lso c alculated s(|w|, w) 



at xq = 0.90, xo = 0.92, xo = 0.95 and have used these results as an input in eq.( ^.17 ). This eventually 
leads to Xq - isolines in the glassy globule phase which are plotted in Fig. 2. 



s(|v|,w) 




3rd virial coefficient, w 



2nd virial coefficient, |v| 



FIG. 1. The dependence of the parameter s (see eq.(5.5)) on the virial coefficients at z = 9 and xo = 



The line, which is associated with eq.(5.22) corresponds to the Ginzburg criterion for fluctuations 
in RS-sector and separates the glassy globule from the liquid one. Obviously, the position of this line 
depends from the value oi £g < < 1 and should be better seen as a crossover from the fluctuating regime 
to the mean field one. In Fig. 2 this line is given at = 0.033. We have not shown more xq - isolines 
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explicitly, but it is important to remind that by changing xq continuously one can span up the whole 
phase digram from left to right. It is interesting that the xq - isolines in Fig. 2 are almost vertical. This 
means that in a real experiment (upon changing the solvent quality |w| by temperature) always some 
particular value of xq is hit in the glassy phase and stays with it as |w| increases. We recall that 1 ~ xq 
is the fraction of replicas which overlap with the strength am = ^m\v'^- 



2.5 



*S 

Is 

o 
u 

-a 

u 

'> 

-a 
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0.5 
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j i 1 x„=0.90 
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FIG. 2. The phase diagram of the polymer globule. The lower solid line corresponds to the coil - liquid globule 
second - order transition, whereas the upper solid line is associated with the first - order liquid - to - solid 
globule transition. Dotted (- dashed) a;o-isolines correspond to glassy state with different values of xo. The other 
a;o-isolines which are not shown here have intermediate values of xq,. 



The fact that on the transition line the value of is less than one, xq < 1 shows that the transition is 
thermodynamically of the first order. This is contrary to p -spin spin glasses [^|-^ and random heteropoly- 
mer [|l^. It is well known that in these cases transition has no latent heat ( i.e. thermodynamically of the 
second order) since xq = \ aX the transition point, whilst the order parameter a undergoes a jump (i.e. 
displays the first - order transition) . In our case the transition is of the first order thermodynamically as 
well as with respect to the order parameter. In Fig. 2 we have shown also the line which corresponds to 
the coil - globule second order transition (see eg. ( 4. 16 )) at the chain length N = 250. The critical value 
|^;|cr on this line is scaled as 1/VlV. 

Fig. 3 shows the reduced density behavior (see eg. ( ^.18 ) at xq ~ 0.88 in the same as in Fig. 2 intervals 
of \v\ and w. It can be seen that the density on the transition line is fairly small to justify the use of the 
virial expansion. 
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c{|v|,w) 




3rd virial coefficient, w 



2nd virial coefficient, |v| ^ 



FIG. 3. The reduced density c as the function of virial coefficients at z = 9 and xo = 0.1 



Let as now calculate the configurational entropy 5conf (or the complexity) which is usually of 
interest for glass forming liquids. The configurational entropy in the liquid globule state can be defined 
as the difference 



Sr.. 



quid 



valley 



(5.27) 



where 



5, 



liquid 



9/i 



RS 



dT 



(5.28) 



and S'vaiicy is the entropy (per particle) which corresponds to one pure state or a valley in the free 
energy landscape. In order to estimate S'vaiiey let us recall from reference ]4^ ] that the order parameter 
a describes the structure of space of valleys through the probability, P{a), that two valleys picked at 
random, have an overlap a. For the 1-RSB - scenario this function has a rather simple form: P{a) = 
XQ6{a) + (1 — XQ)6{a — Um)- The dimension of the matrix Aq^, which within 1-RSB is parameterized by 
n(x{) — l)/2 parameters a, becomes negative at n ^ and <_^o 5: 1- Therefore the free energy /rsb 
(see eq.(p!^)) becomes negative, and the factor — (1 — xq) in eq.( |3.4[ ) can be treated as "fraction of valleys" 
with overlap ■ I n th is case the intra- valley free energy can be obtain by dividing out a common factor 
of —(1 — xa) in eq.(5^): /vaiicy — — /rsb/(1 — a;o). Then we can estimate S'vaiiey in a following way: 



1 



valley 



RSB 



1 - Xo dT 



(5.29) 



The configurational entropy in the glassy globule state is given by the same way as in the p -spin spin 
glasses Hpf: 



•Sconf = [l/'(l) - ■0(1 - 2;o)] , 



(5.30) 



where ^{x) is the digamma function. 
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FIG. 4. The configurational entropy 5conf(r) as the function of temperature whilst the globule is cooling from 
the liquid to the glassy state 



conf undergoes a 
simulations. It 



By cooling the system alo ng th e x g - is oline at xq = 0.88 , we have calculated the configurational 
entropy for hquid (see eqs. ( 5.27 ) - ( 5.29 )) and glassy ( eq.( 5.30| ) globule. Figure 4 gives the result 
of this calculations. As discussed above the transition is of the first order, i.e., the S*, 
jump which is qualitatively in line with Molecular dynamics |^,^ and Monte Carlo |2^ 
must bear in mind that the transition temperature in Fig. 4 is also defined by the Ginzburg criterion for 
fluctuations (see the intersection point of the upper solid line and the xq — 0.88 isoline in Fig. 2), so that 
this transition can be treated as a crossover from the fluctuating regime to mean field, where only few 
states are dominated. This possibility is shown in Fig. 4 by dotted line. On the other hand, this behavior 
is quite different from low molecular systems or polymer melts where Sca-di goes to zero continuously as 
soon as T ^ Tk- 



VI. CONCLUSION AND OUTLOOK 



We have shown that the homopolymer globule problem can be formulated within the n — > limit of 
free n-component V'f j '06 " field theory and then be mapped onto the disordered one - component model 
with a non - Gaussian random noise. It is of interest, that the statistical moments of this noise can be 
expressed only through the virial coefRcients of the pure model, the fact which gives grounds to discuss 
the self - generated disorder. On the other side, the n limit is responsible at the same time both for 
the polymer conformations and for the nontrivial structure of the replica space. Physically this means 
that units with competitive interactions and stipulated by the connectivity posses good pre-conditions 
for the glass formation. 

The Legendre transformation from the pair fields ipaif)''l'bif) to the Parisi overlap parameters Qabif) 
leads to the effective Hamiltonian which is akin to corresponding expressions for the spin glass models. 
We have investigated this resulting replicated model within mean field Landau - type treatment. First of 
all, the RS - solution of the corresponding equation is associated with the conventional coil - liquid globule 
transition whereas the RSB - solution in the deeper globule state is related with the glassy regime. This 
mean field glassy globule phase can be only assured if the fluctuations in RS - and RSB - sectors are small. 
In this paper we have only studied the fluctuations in the RS - sector and have sketched the corresponding 
Ginzburg criterion line on the phase diagram. We have calculated the configurational entropy for the 
liquid and glassy globules and shown that the transition between them is a first order one, as also is in 
the MD and MC computer simulations |^,^,^ , whereas for the random heteropolymer the freezing 
is thermodynamically a second order transition. On the other side, the computer simulation [ p3| and 
experiment [p4[ shows that the folding in proteins has a latent heat, i.e. the random heteropolymer is a 
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poor model of protein - folding thermodynamics. It was shown that in the heteropolymer models 

where some sequences of monomers have especially law energy in their native conformation the folding 
is thermodynamically a first - order transition. 

Nevertheless, our calculations shows that already for the homopolymer globule, because of competitive 
interactions (virial coefficients of different signs) and connectivity , the transition with a latent heat is 
possible. In this respect the revisiting of the random heteropolymer problem within our field theoretical 
approach (as oppose to the density functional method [T^-[l^) would be very interesting. This could 
elucidate the problem how the self - generated disorder interplay with quenched disorder and eventually 
modify the freezing conditions. 

The dynamical aspects of the thermodynamic transition discussed above is a matter of crucial impor- 
tance. The first and foremost question which should be investigated is related with the formation of the 
entropic droplets [^p^] in the RSB - sector of the replica space. 
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APPENDIX A: FIELD-THEORETICAL REPRESENTATION FOR A GENERAL 

SELF-INTERACTING CHAIN 

Let us represent the Hamiltonian of the self-interacting chain in the following form 

N 2 


N 

oo „ m 

+ E T;^^ / dsodsi ...ds^Y[S (r(so) - r{sij) , (Al) 

m=l ^ Q i=l 

where Vm+i denote the virial coefficients. It is convenient to introduce the density 

N 

p(r)= / dsS{r-ris)), (A2) 







then the partition function of the polymer chain 

Z = J Dr{s)exp{-H[r{s)]} (A3) 



takes the form 



N 



Z = j Dr{s)Dp{r)5 p{v) ~ j dsS (r - r(s)) 



X exp ■ 



N „ 
d , Vm+l 



2a2 y \ds ^ (m + 1)! 

™=i 



By making use of the integral representation for the (5 - function cq.(A4) can be recast in a form 
Z = y" DY{s)D(\){v)Dp{r)&^\>\^ J d'^r(j){r)p{r) 

- E ^^^^ J rf'^P^+'wj / d''rdVG(r,r';[0];iV), (A5) 



17 



where 



G(r,r';M;7V) - j i?r(s)exp J - A j ds j dscjy{v{s)) 



(A6) 



r(0)=r 

The corresponding equation for G reads 



G (r,r'; [</)]; TV) =,5(r-r')<5(iV) 
It is convenient to make the Laplace transformation 



G(r,r';[0];A*) = j dTVG (r, r'; [0]; iV) exp(-/.7V), 



after which the equation for G (r, r'; [0]; /i) yields 

/i:>VV'(r)V'(r')exp|-i/d'*rV'(r) M " §3' 



(A7) 



(A8) 



G(r,r';[<^];/x) = 



(/.(r)] ^(r')} 



(A9) 



In order to avoid the denominator in eq. (|A9|) we should upgrade the field ■;/; by introducing n-component 
field Tp = {'01) '02j • ■ • V'n}- Then by using the replica trick, we have 



G (r, r'; [0]; ^) = Urn y \{ 7^^„Vi(r)0i(r') exp | - ^ y d'*r^a(r) 



Ai- — V2 + z0(r) 



0a(r') (AlO) 



One can make now the Laplace transformation from both sides of eq.(A5) and substitute in it eq.( A10| ). 
The integration first over the field 0(r) and then over p(r) results in the following expression for grand 
canonical partition function of a polymer chain with ends fixed at points r and r': 



/ ■a 
Jl £>i/'QVi(i")'0i(r') exp |-i?„ 



where 



0;^ 



"'if 



-1 a— 1 



m+l 



(All) 



(A12) 



APPENDIX B: THE CONNECTION TO A RANDOM MODEL 



Here we prove that in a poor solvent the polymer problem (see eqs.(2.3) and (2.4)) can be connected 
to a one- component random model. We shell restrict our consideration only to second and third virial 
coefficients. The generalization to the case with an arbitrary number of virial coefficients is straightfor- 
ward. 

Let us consider the one component random model with Hamiltonian 



(Bl) 



where t{r) is a non-Gaussian random field with the generating functional 
${p(r)} = J Dt{r)P{t{r)} exp 1^- J d'^rt{r)p{r 



exp 



d'^rp^{r) - 



3!5 



d'*rp3(r) 



(B2) 
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In eq.(p32|) P{t{r)} is the distribution functional of the field t{r). One can easily check that the replication 
of the Hamiltonian (Bl) and the subsequent averaging over t{r), i.e. 



/n ( n ^ 

a=l y 0=1 J 



(B3) 



leads to the effective replicated Hamiltonian (2.4). 

^From the probabilistic interpretation (B2) one can explicitly find the central moments of t{r). The 
expansion of both sides of eq.(B2) yields 

(-1)' 



E 

oc 



m! 

k 



dri... drmt{ri)t{r2) ■ ■ ■ t{r,n)p{ri)p{r2) ■ ■ ■ p{r„i) 



fc!(-l) 



k=o 1=0 ^ ' ^ 



W 

3!8 



(B4) 



By making the m-terms of p{r) on both sides of eq.(B4) equal, one gets 



/ 



dri... drrat{ri)t{r2) ■ ■ ■ t{r,n)p{ri)p{r2) ■ ■ ■ p{rm) 



[m/2] 

Em! 
(m- 2fc)!(3fc- m)! 

fc={m/3} 



/" 3k~~-Tn r /t 
d'rp\r)\ J d'^rpHr) 



'2k 



(B5) 



where [m/2] stands for the greatest integer number less then m/2 and {m/3} is the smallest integer 
number larger then m/3. By making use the representations Jd'^rp^(r)= J d'^rid'^r2S{ri — r2)p{ri)p{r2) 
and J d'^rp'^{r) = J d'^r3d'^rid'^r^5{Y^ — V4)5{r^ — V3)p[vT,)p{vi)p[r^) one finally gets: 



V2] 



t{T^)t{T2)...t{r^)= 



(m- 2fc)!(3fc- m)! V 8 



Zk — Trl 



W 

3!8 



%-2k 



fe={m/3} 

X (5(ri - r2)^(r3 - r4) . . . 5{YQk-2m-l - r6fe-2m) 

2m+l I'e/c — 2m+2 2m+l — 2m+3j ■ ■ • ^\^7n — 2 I'm— 1 ^m — 2 I'm 

). (B6) 



In eq. (p36| ) the second line includes 3k — m i5-functions with arguments successively pairwise divided 
between 6fc — 2m points. The third line includes m — 2k (5-functions so that in each successive pairs of 
them one r - point is common. 

Let us consider some particular cases of eq.(B6). 

i)m = 1. Then {m/3} = 1, [m/2] = and 

= (B7) 
n)m = 2. Then {m/3} = 1, [m/2] = 1 and eq.@ reads 

U,l 

(B8) 
(B9) 
(BIO) 
(Bll) 



t{ri)tir2) ^ ^S{r, ^ r2) 

i) At m — 3, {m/3} — 1, [m/2] — 1 (i.e. k=l) and one gets 

w 



t{ri)t{r2)t{r3) = -(5(ri - r2)(5(ri - rg) 



iv) At m = 4, {m/3} = 2, [m/2] = 2 (i.e. k=2) and 

3 



t(ri)t(r2)i(r3)t(r4) = -^H'^Sin ~ r2)(5(r3 - 
16 



v) Finally at m = 5, {m/3} = 2, [m/2] = 2 and 

5! , 



t(ri)t(r2)t(r3)t(r4)t(r5) 



3!64' 



u6{ri - r2)^(r3 - r4)(5(r3 - rs) 



The important feature of these moments is that all of them are positive, which means that t(r) is real. 
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APPENDIX C: SPATIAL FLUCTUATIONS FOR WEAKLY INHOMOGENEOUS GLOBULE 



In this appendix we give the Landau-expansion only for the case when q and coefficients T^^) , r^^) , r^'') 
are weakly k - dependent in the RS - sector. The spatial Fourier transformation in eq.(3^) leads to the 
following effective Hamiltonian 



HRs{qm 



where 



^r(2)(k)g(k)g(-k) 



{2n 

(27r)3d 



d'^kid'^k2 



r(3)(ki,k2)q(ki)g(k2)g(-ki -k2) 



(27r)2d 

r(4)(ki,k2, k3)g(ki)q(k2)g(k3)<z(-ki - k2 - ks) 



d'^K 1 

[i-^+M] [i(.^-k)2+^] ' 



(Cl) 



(C2) 



r(^^(ki,k2) 



w 



d'^K 1 

(2^ [i'^' + [U^ - ki)^ + - ki - k2)2 + 



(C3) 



rW(ki,k2,k3) 



For weak fluctuations around MF - solution to be investigated one should estimate Hessian 



5q{ri)5q{Y2) 



-2r(2)(ri-r2) + 6g„, / d^r^T^^Xv^ - r^,v^ - v^) 

q=qm 

+ 12ql J dSrf''r4rW(r2-ri,r3-ri,r4-ri) + ... , 



or in the Fourier space 
6q{ri)6q{r2) 



= 2r(2) (k) + 6<z,„r(3) (k, k = 0) + i2g2^r(4) (k, k = o, k = o) 



(C4) 



(C5) 



(C6) 



The resulting effective Hamiltonian expansion around the MF - solution takes the following form 



1 f d I 
ffRs{q™,A(7(k)} = ifR,s{g™} + - J 



I f d-'k S^Hns 



\'i 5q{'k)dq{-k) 



Ag(k)Ag(-k) 



9=9,7 



= Aql-^ + Bql + Cqi 
d'^k 



xAg(k)Aq(-k) 



r(2) (k) + 3g„r(3) (k, k = 0) + 6g^r(4) (k, k = 0, k = 0) 



where 



and 



Ag(k) = g(k) - q„ 



A = r(2)(o) , B = r(3)(o,o) , c = r(4) (0,0,0) 



(C7) 
(C8) 
(C9) 



Th e ex pressions for A, B and C are given by eqs. ( p.7[ ) -(^^). It is easy to estimate the integrand in 
eq.(C7) at small k (weak inhomogeneity) . The straightforward calculations yields 
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(2^ 



{Xi\vlw) + ikafY{\vlw)}Aq{k)Aq{-k), (CIO) 



where 

X = A + 3g„,B + 6g^C (Cll) 

and 

~ [2TTda^ ) Uf,^ \ 2 j Iv27rda2 j ' ^ 



9^- 



From (CIO) it is obvious that the corresponding correlation function 



lAg(k)l^) ^ ^ . (C13) 

/ X{\v\,w) + {ka)'^Y [\v\,w) 
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